Ring-like spin segregation of binary mixtures in a high-velocity rotating drum 
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We present molecular dynamics simulations on the segregation of binary mixtures in a high- 
velocity rotating drum. Depending on the ratio of particle's radius and density, the ring-like spin 
segregation patterns in radial direction show similarities to the Brazil-nut effect and its reverse 
form. The smaller and heavier particles are easier to accumulate towards the drum wall and the 
bigger and lighter ones towards the drum center. Furthermore, we quantify the effects of particle's 
radius and density on the segregation states and plot the phase diagram of segregation in the pb / ps 
- r^jra space. The observed phenomena can also be well explained by the combined actions, i.e., 
the percolation effect and the buoyancy effect. 
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I. INTRODUCTION 

Granular materials are widespread in nature and daily 
life such as rocks, sands and sugars. The research on such 
a group of discrete particles is a particularly intriguing 
subject for engineers and physicists because of its gas, 
liquid and solid like properties [l), [1] . A common exam- 
ple of segregation of granular mixtures usually happens 
ranging from industrial applications such as pharmaceu- 
tical processing and mineral transportation to natural 
phenomena such as snow avalanche and sand sedimenta- 
tion in riverbeds [3i, i^]. In the last decades, there have 
been a number of researches to control and characterize 
this peculiar property of granular mixtures. It is found 
that granular mixtures can segregate under a variety of 
conditions, i.e., shear [H, vibration [6], and rotation Q. 
The results show that the segregation states are governed 
by particle size [1] j_density Q , friction coefficient [l0| , vi- 
bration intensity [ll| and even the pressure gradient of 
air in the interstices between particles [l^l ■ Though large 
progresses have been made on this phenomenon, the be- 
hind physics remains still interesting and new behaviors 
or effects are to be found and investigated. 

In general, the spontaneous segregation always hap- 
pens during the granular flowing process. Thus the sys- 
tem of granule flowing in a horizontal rotating drum is 
a good candidate for the experimental and theoretical 
convenience Among these studies, two special 

cases of segregation of binary mixtures, where the parti- 
cles differ in size (known as S-type system) and density 
(known as D-type system), have received intensive atten- 
tions [l^ [l^ . In an S-type system, the smaller particles 
are more easily sieved down though the interstices be- 
tween particles, which is always called as the percolation 
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effect. On the other hand, in a D-type system, the par- 
ticles with higher density firstly sink down contrasted 
to the lighter ones, which is always named as the buoy- 
ancy effect. Furthermore, due to the combined actions of 
percolation and buoyancy effects in same direction, the 
degree of segregation can be reinforced in an S+D-type 
system, in which the bigger particles have a lower density 
and the smaller ones have a higher density. Reversely, the 
segregation state can be changed completely in an S-D- 
type system, in which the bigger particles have a higher 
density and the density driven effect takes an absolute 
role. 

Many reported studies found that the rotational ve- 
locity is also a key control parameter on the segregation 
pattern [13, [3| • Depending on the rotational velocity of 
the drum, different flow regimes can be obtained, such 
as avalanchin g, r olling, cascading, cataracting, and cen- 
trifuging [Tol, 120|. Especially, Nityanand et al. found 
that the smaller particles accumulate in a classical core 
of the bed for the rolling regime [2l|. Such segregation 
can disappear when the granules flow in the cascading 
regime. Then the reversed segregation begins to occur 
when the granules flow in the cataracting regime. The 
segregation will become more pronounced under the cen- 
trifuging regime, in which the smaller granules adjacent 
to the drum wall are well graded and those having larger 
size far from the wall flow coarsely. To our knowledge, no 
much attentions associated with the centrifuging segre- 
gation of binary mixtures have been carried out. A more 
detailed study on the centrifuging granular flow in a ro- 
tational drum would be useful for a better understanding 
of the segregation of binary mixtures. 

In this paper, we make a simulational study of the seg- 
regation of particles with different sizes and densities in 
a high-velocity rotating drum. In Sec. II, the simulation 
model is described and the segregation index is intro- 
duced to deflne the degree of segregation. In Sec. Ill, 
a series of systems, i.e., S-type, D-type, S-t-D-type and 
S-D-type systems, are simulated and the different ring- 
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like spin segregation patterns are given. In Sec. IV, we 
quantify the effects of particle's radius and density on 
the segregation states and plot the phase diagram in the 
Phi Ps - Thl"fs space. Then a simple explanation consistent 
with the findings is discussed. Section V summarizes the 
conclusions of the work. 



II. SIMULATION MODEL 

The simulations are performed in a quasi two- 
dimensional drum as shown in Fig|T] The drum is 
D = 200.0 mm in diameter and St — 1.0 mm in thickness. 
In this paper, all runs are begun from a randomly mixed 
state. The particles are randomly created according to 
their relative volumes in the upper part of drum and fall 
under the self-gravity when the rotational velocity is set 
to zero. The gravitational acceleration in the radial plane 
of drum is set to g = 9.8 m/s^. Equal volumes of the two 
types of particles are used in all simulations. Depending 
on the rotational velocity, different flow regimes can be 
produced. Two indicated rotational frequency, fn and 
Jl , are critical for marking the stable fiow bands of par- 
ticles [12, [2^ . If the rotational frequency is much higher 
than fu, all particles form a solid layer adhering to the 
drum wall due to the larger centrifugal force and finally 
rotate together with the drum. In this case, the par- 
ticles have no relative motion and no clear segregation 
appears. If the rotational frequency is less than f^, the 
majority of particles can not obtain enough energy to 
form a centrifuge flow. Here, /h and /l are values being 
about 100.0 rmp and 200.0 rmp, which are also found 
to be closely correlated with the filling fraction and the 
friction coefficient between particles and between parti- 
cles and drum wall. When the frequency is lowered to a 
moderate value between fn and /l, some particles con- 
tinue flowing along the drum wall, while other particles 
detach from the wall and flight as a free falling body. 
In the simulations, the filling height of particles is set 
to Hf = 0.6D. The drum rotates clockwise around its 
center at a moderate high frequency, / = 150.0 rpm. 
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particle, which has been used in our previous studies 
[24l [25! . In a simulation time step, the positions and ve- 
locities of particles are updated by integrating the New- 
ton's second law. To form a quasi two-dimensional sys- 
tem, the particles are circular in the radial direction of 
drum and the particle's thickness in the axial direction 
of drum is of the same value of drum thickness. In the 
simulation, both the translational motion in the radial 
plane of drum and the rotational motion in the axial di- 
rection of drum are considered. The force between two 
contacting particles is calculated in normal and tangen- 
tial directions. The Kuwabara-Kono model is used to 
describe the normal interaction 129, [23 , 



(1) 



The tangential component is determined by the minor 
of the tangential force with a memory effect and the force 
of dynamic friction. 



F[. = -mm{ktQj,p.Fl 



(2) 



In Eqs.(IT]) and ([2]), i and j are the indexes of the par- 
ticles. S^ij — max(0,(iy — jx^ — Xj|) is the overlap of 
the contacting particles i and j, d is the sum of radiuses 
of the contacting particles and rj. Cij denotes the 
displacement in the tangential direction that took place 
since time to, when the contact was first established, i.e., 
Qj{t) = Vs{t')dt'. Vij = V'-Bn + V^jBr is the relative 
velocities of two particles. The indexes n and r represent 
the normal and tangential directions at the contact point, 
respectively. kt and ?7„ characterize the stiffness and 
damping of the granular materials and are related to the 
collision time t„ and the masses of contacting particles rrii 
and rrij. The detailed values can be calculated as follows. 
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E is the Young modulus, v 



the Poisson ratio. The coefficient of sliding friction is set 
to /i = 0.5. The collisions between the particles and the 
wall of drum are treated as particle-particle collisions, ex- 
cept that the drum has an infinite mass and a diameter 
D. 

In the simulation, the center position of the masses is 



— J2i=i ^i- The angular velocity 
'^^Ji uji. Here, p is the species 



FIG. 1: A snapshot of system at the beginning of simulation. 



easily calculated, hp 

of the masses is ujp 
of particles and Np is the total number of p species of 
particles. Zi is the distance of particle i with respect to 
the edge of drum. uJi is the angular velocity of particle i. 
In order to have an quantified description of the degree of 
segregation, the segregation index x is usually introduced 



Molecular dynamics simulation is adopted by using the 
discrete element method to describe the motion of each 



X = 2.0 



hb + hs 



(3) 
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In Eq. ([3]), /ib and hs are the center positions of the 
bigger particles and the smaller particles, respectively. 
As a comparison with the segregation in a vertically vi- 
brated container , we also call the segregation state 
as the Brazil-nut effect (BNE) , in which the bigger parti- 
cles migrate towards the center of drum and the smaller 
ones towards the drum wall. If the segregated state is re- 
versed, i.e., the bigger particles accumulate towards the 
drum wall and the smaller ones towards the center of 
drum, the corresponding segregated state is so named 
the Reverse Brazil-nut effect (RBNE). Apparently, the 
BNE happens for a positive while the RBNE for a 
negative x according to the Eq. 



III. RING-LIKE SPIN SEGREGATION 
PATTERN 



The simulations are firstly performed for an S-type sys- 
tem where the particles are of the same density but differ 
in size. Figl^Ja) shows the segregation pattern of the bi- 
nary mixtures after 500 rotations for rb = 1.5 mm bigger 
particles (red circles) and = 1.0 mm smaller particles 
(blue circles). We can see that two ring-like spin seg- 
regated patterns in the radial direction appear except a 
mixed region near the drum wall. The bigger particles 
migrate towards to the center of the drum and form a 
ring-like layer with a smaller radius, while the smaller 
particles surround at the periphery. In addition, the in- 
ner ring with bigger particles is actually a semi-circle, 
in which a few of bigger particles leave from the flowing 
layer and flight freely at the upper right region. 

To give more insight on the segregation process of dif- 
ferent species of particles, we plot the evolutions of the 
center position and the angular velocity of bigger parti- 
cles and smaller particles as shown in Figl^Kb). Figl2](b) 
shows that the segregation happens in the first hundred of 
rotations, i.e., the bigger particles are pushed towards the 
center of drum {hb — 22.06 mm) and the smaller particles 
towards the drum wall (hg = 15.12 mm). At the mean- 
time, the averaged angular velocity of both the bigger 
particles and the smaller particles increases very fast and 
approaches nearly to that of the drum, ujb = 15.45 s^^ 
and ujs = 15.67 s~^. Further calculation shows that the 
particles at the mixed region just adjacent to the drum 
wall have the same angular velocity as that of the drum. 
The particles have no relative motions and this results in 
that the particles in the mixed region will not segregate 
even if the rotation continues. 

Second, the simulations are carried out for a D-type 
system in which the particles are of the same size (r^ = 
r; = 1.5 mm), but the densities of heavy particles and 
light particles are set to ph = 3.0 x 10'^ kg/m^ and 
pi — 1.0 X 10^ kg/m^, respectively. The particles also 
segregate into two ring-like spin layers by density as 
shown in Fig|31[a). The two rings have high concen- 
trations with heavier particles locating at the periphery 
and lighter ones at the inner side. Likewise, the inner 
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FIG. 2: (color online) (a) Image of the S-type system after 
500 rotations. The densities of both bigger particles (open 
symbols) and smaller particles (solid symbols) are of the same 
value. Pi, = ps = l.Ox 10^ kg/m^. The radiuses of particles are 
for different sizes being Vb = 1.5 mm and Vs = 1.0 mm, respec- 
tively, (b) Evolutions of the center positions (solid symbols) 
and the angular velocity (open symbols) of bigger particles 
and smaller particles corresponding to the system of (a). The 
square symbols and circle symbols are results for the bigger 
particles and smaller particles, respectively. 



ring is a semi-circle, where some lighter particles detach 
from the flowing layer on the upper right region. As 
shown in FiglH^b), the segregation also proceeds quickly. 
The heavier particles migrate towards the drum wall 
{hh — 11.24 mm) while the lighter ones towards the drum 
center [hi — 31.89 mm). Furthermore, the heavier par- 
ticles approaching closer to the drum wall have a higher 
angular velocity (cj/i — 15.58 s~^) compared to the lighter 
particles (w/ = 12.99 s^^). 
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FIG. 3: (color online) (a) Image of the D-type system after 
500 rotations. The radiuses of both heavier particles (solid 
symbols) and lighter particles (open symbols) are of the same 
size, rh = ri — 1.5 mm. The densities of particles are for 
different values being ph — 3.0 x 10^ kg/m^ and pi = 1.0 x 
10^ kg/m^, respectively, (b) Evolutions of the center positions 
(solid symbols) and the angular velocity (open symbols) of 
bigger particles and smaller particles corresponding to the 
system of (a). The circle symbols and square symbols are 
results for heavier (lighter) particles. 

A more complicated segregation happens for the cases 
where the binary mixtures are different in both size and 
density. We set fixedly the radiuses of bigger particles 
rb — 1.5 mm and smaller particles — 1.0 mm, and then 
change the density ratio of bigger particles to smaller 
ones. FiglDJa) shows the segregation pattern after 500 
rotations for an S+D-type system, where the densities of 
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bigger particles and smaller particles are the values being 
Pb = 1.0 X 10^ kg/m^ and Ps = 3.0 x 10^ kg/m^, respec- 
tively. Compared to the segregation pattern in FiglJJa), 
we can see that two more clearer ring-like spin segrega- 
tion layers appear, in which the smaller and heavier par- 
ticles gather at the periphery and the bigger and lighter 
particles at the inner side. If the density of bigger parti- 
cles is increased larger than that of smaller ones to form 
an S-D-type system, the segregation can be reversed com- 
pletely, i.e., the bigger particles appear at the outer side 
and the smaller ones at the inner side. FigUJ^b) shows a 
typical RBNE segregation pattern, where the density of 
bigger particles is set to pb — 3.0 x 10'^ kg/m^, and that 
of smaller ones = 1.0 x 10"^ kg/m"^. 




FIG. 4: (color online) Combined size and density segrega- 
tion patterns of binary mixtures after 500 rotations. The ra- 
diuses of bigger particles (open symbols) and smaller ones 
(solid symbols) are for different sizes being rt — 1.5 mm 
and Ts = 1.0 mm, respectively. (a) An S+D-type sys- 
tem. The densities of bigger and smaller particles are set 
to Pi = 1.0 X 10^ kg/m^ and ps = 3.0 x 10^ kg/m^, re- 
spectively, (b) An S-D-type system. The densities of bigger 
and smaller particles are set to pf, = 3.0 x 10^ kg/m'' and 
ps = 1.0 X 10"^ kg/m^, respectively. 



IV. DISCUSSIONS 

To give a quantitative description of the above ring- like 
spin segregation pattern, we first plot the dependence of 
the segregation index on the density ratio of Pb/ps for 
different values of radius ratio rb/rg- In FiglH at a given 
radius ratio rt/rs, the system changes from the S-|-D-type 
to the S-D-type when the density ratio of Pb/ps grows. 
The segregation pattern of binary mixtures moves from 
the BNE to the RBNE. Furthermore, as shown in FiglSJ 
the whole curve shifts upwards as the radius ratio rb/rg 
increases, which means that the switch of the BNE to 
RBNE needs a higher value of density ratio of Pb/ps- 

Figl6] summarizes the above findings in the phase dia- 
gram of Pb/ Ps-fb I space by using the segregation index. 
As shown in Figll the BNE and the RBNE are indi- 
cated by the symbols of the upward-pointing triangles 
and downward-pointing triangles, respectively. The solid 
lines give the boundaries of the BNE and the RBNE, 
where the absolute values of x ^-re actually very small 
and the corresponding states are almost mixed. In Fig|51 
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FIG. 5: Dependence of the segregation index, X) on the ra- 
dius ratio of pb/pa- The radiuses ratio of r^/rs are values 
being 1.5 (square symbols) and 2.0 (circle symbols), respec- 
tively. As pb/ Ps increases, the system moves from a BNE to 
a RBNE segregation pattern. The dashed line corresponds to 
the boundary of BNE and RBNE. The solid line is guided for 
eyes. 



it can be easily seen that the BNE appears at the upper 
left region and the RBNE occurs at the lower right re- 
gion. Furthermore, the boundary shifts right when the 
density ratio of bigger particles to smaller ones increases. 
This result suggests that the BNE is favored for the lower 
Pb I Ps and the larger rb/rg, while the RBNE for the higher 
Pb/ Ps and the smaller Vb/rg- 
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FIG. 6: Phase diagram of the binary mixtures in the pb/ps 
- Tb/vs space. The upward-pointing triangles and downward- 
pointing triangles represent the segregation patterns being the 
Brazil-nut effect and the Reverse Brazil-nut effect. The solid 
lines show the boundaries of two segregation patterns. 

The above segregation phenomena due to size and den- 
sity might seem to result from the two effects, i.e, per- 
colation effect and buoyancy effect. At the case of the 
centrifugal flow regime, the centrifugal force is always 
in the radial direction and the outermost region, i.e., the 
drum wall, is the place of bottom region of granular flow- 
ing layer. For an S-type system, the smaller particles are 
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easier to be captured by the voids between the particles 
in the flowing process and migrate to the bottom region. 
On the other hand, for a D-type system, the centrifu- 
gal force is linear to the mass of particle, i.e., oc rup. 
It means that the particle of higher density possesses a 
larger centrifugal force compared to that of lower den- 
sity. To minimize this centrifugal potential of system, 
the species of particles with higher density are more likely 
pushed to the drum wall, i.e., the bottom region. This 
kind of segregation effect due to the difference of den- 
sity is here analogous to the buoyancy effect happened 
in the vertical vibrated binary mixtures. According to 
the discussions of percolation effect and buoyancy effect, 
it is reasonable to see what happened in the S+D-type 
system and the S-D-type system. For an S+D-type sys- 
tem, the percolation effect and buoyancy effect act in 
the same direction and this combined effects result in a 
more segregated state. In reverse, the percolation effect 
and buoyancy effect will compete with each other in an 
S-D-typc system. Thus the final state of binary mix- 
tures might be segregated or mixed, which depends on 
the combined actions of the percolation effect and buoy- 
ancy effect. For example, when the percolation effect 
overcomes the buoyancy effect, the BNE happens, while 
when the latter dominates, the RBNE occurs. If the two 
effects are almost equally, a mixed state will finally yield. 

V. CONCLUSIONS 

In this paper, a ring-like spin segregation pattern of bi- 
nary mixtures in a high- velocity rotational drum is inves- 
tigated by computer simulations. The radial segregation 



we observed here is analogous to that a vertical vibrated 
container, i.e., the Brazil-nut effect, in which the larger 
particles accumulate towards the center of drum and the 
smaller ones towards the drum wall, and the Reverse 
Brazil-nut effect, in which the larger particles migrate 
towards the drum wall and the smaller ones towards the 
center of drum. The dependence of the segregation index 
on the density ratio of Pb/ Ps further quantify the switch 
between the BNE and The RBNE. The phase diagram in 
the Pb/ Ps - rb/vs space shows that the BNE is favored for 
lower density ratio and larger size ratio of bigger particles 
to smaller ones and the RBNE for higher density ratio 
and smaller size ratio of bigger particles to smaller ones. 
We also propose a simple analysis that the simulated phe- 
nomena result from the combined actions of percolation 
effect and buoyancy effect. When both effects act in same 
direction, the BNE happens, while when the density ra- 
tio of bigger particles to smaller ones is increased high 
enough and the buoyancy effect exceeds the percolation 
effect, the RBNE occurs. Since the segregation is of very 
importance in processing granular mixtures, the results 
in this paper have potential industrial and theoretical 
values. 
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